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ABSTRACT 

This paper presents a study of the dynamic behavior of a rotor system supported by 
auxiliary bearings. The steady-state behavior of a simulation model based upon a 
production jet engine is explored over a wide range of operating conditions for varying 
rotor imbalance, support stiffness and damping. Interesting dynamical phenomena, 
such as chaos, subharmonic responses, and double-valued responses, are presented and 
discussed. 
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NOMENCLATURE 

Cq — auxiliary bearing support damping, lb. s/in. 

F n = normal force, lb 

F X = external force vector acting on the rotor in X direction 
Fy = external force vector acting on the rotor in Y direction 
I a = rotor inertia matrix 

Kb — auxiliary bearing support stiffness, lb/in. 

Kq = contact stiffness, Ib/in. 

Mb = auxiliary bearing mass, lb.s 2 /in. 

Mk — mass of kth rotor element, lb.s“/in. 

N = total number of modes considered 
Nsi = node number at auxiliary bearing #1 
N B2 — node number at auxiliary bearing #2 
Qx = rotor modal coordinate vector in X direction 
Q Y = rotor modal coordinate vector in Y direction 
Rq = radius of auxiliary bearing bore, in. 

Rr = radius of rotor journal, in. 

Xh = rotor physical coordinate vector in X direction 
Y fi = rotor physical coordinate vector in Y direction 
e = rotor imbalance eccentricity, in. 
g = gravitational acceleration, in./s 2 
t — time, s 

A = deformation at the contact point, in. 

r = ^ T i a ^ 

= rotor free-free modal rotation matrix 
f) = rotor operating speed, rad/s 

= rotor free-free modal displacement matrix 
6 = Rq — 7? /?, auxiliary bearing clearance, in. 

■aECfiMNCi FAGE BLANK NOT BLMEb 
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( = modal damping coefficient 


INTRODUCTION 

Active magnetic bearings are one of the most innovative recent developments in the turbomachinery 
field. This technology provides the potential for significant improvements over other types of rotor 
support, including elimination of wear and bearing friction-related energy losses as well as a means 
of actively suppressing rotor vibration. A critical component of any magnetic bearing design is the 
auxiliary bearing, which protects the soft iron core of the magnetic bearing and provides rotor support 
in case of overload or failure of the primary (magnetic) bearing. Magnetic bearing systems appear 
to provide particularly great promise for use in aeronautical applications. In this regard, current 
effort is directed toward developing jet engines with rotors supported by magnetic bearings. For such 
applications, safety is an important concern. Toward this end, it is desirable to design the system to 
operate with auxiliary bearing support for an extended period of time. 

A number of different bearing types have been suggested as auxiliary bearings. These include bushings, 
rolling element bearings, and various types of journal bearings. The most commonly considered are 
rolling element bearings. The major disadvantage associated with using rolling element bearings (or 
bushings) is the requirement of a clearance between the rotor and the inner race of the bearing, without 
which many of the advantages associated with using magnetic bearings would be reduced or eliminated. 
This clearance introduces a nonlinear dynamical feature which may significantly impact the behavior 
of the rotor. 

There are quite a number of studies in the literature concerned with nonlinear rotordvnamics. 
Yamamoto (1954) conducted a systematic study of rotor responses involving bearing clearance effects. 
Black (1968) studied the rotor/stator interaction with a clearance. Ehrich (1966, 1988 and 1991), Bently 
(1974), Muszynska (1984) and Childs (1979 and 1982) observed and studied subharmonic responses 
associated with clearance effects. While this work has served to greatly enhance the understanding 
of such systems, more detailed study is needed. Much of this earlier work was conducted from the 
perspective that the clearance is a result of manufacturing error or misfitting and is best eliminated. 
However, in a rotor system fitted with magnetic bearings and auxiliary bearings, the clearance becomes 
a design parameter rather than an irregularity. From this point of view, it is important to develop a 
better understanding of the expected dynamic responses. Such knowledge will provide guidelines for 
the selection of auxiliary bearing parameters. 

There has been relatively little work available in the open literature that is specifically concerned with 
auxiliary bearings. Two papers that are directly related to research on auxiliary bearings are Gelin 
et al., (1990) and Ishii and Kirk (1991). Both of these papers are mainly concerned with transient 
responses. The current work is specifically concerned with developing a better understanding of the 
expected steady-state dynamical behavior for an auxiliary bearing supported rotor. Simulation results 
for a rotor based upon a production jet engine are presented and discussed. 

The model used in the current study has two principal components - the rotor and the auxiliary 
bearings. The rotor is modelled using the free-free bending mode shapes and natural fiequencies 
obtained through finite element analysis. The finite element code uses 31 stations and the first four 
modes (two rigid body and two flexible modes) are included in the simulation model. Figure 1 shows 
a schematic diagram of the finite element model. The rotor equations of motion can be expressed in 
terms of modal coordinates as: 


Qx+'-C^nQjY + ftTQv + u 'nQ.V 

+ 2Sl(u n Q y = ^ T F.v , U- a ) 

Qy+^UJnQy — OTQ.v t V 

- ’ift^WnQ.Y = <& r F y. (Fb) 



398 




FIGURE I DIAGRAM OF THE FINITE FIGURE 2 AUXILIARY BEARING MODEL 

ELEMENT MODEL 


where 


F.v = {F. Yl> F\2, ■ • 

■yF Xm }- T 

Fy = {Fy\, Fy2, ■ ■ 
Qy = * _1 Xr, 

Qy = 

.,Fy m }- T , 


/w ni 0 

... 0 x 

u; n = 

0 w’ nj 

... 0 


Vo o 

• - - w n n ‘ 

Xfi = {A'fij, A R 2> • 

} _T 


Y« = {Vm.i'/w r Rm }- r . 

(m = total number of nodes) 

The physical displacements of the rotor at the two auxiliary bearing locations can be obtained using 
the following coordinate transformation: 


;V 

■^Rk = ^ Qxiy 
»=1 

.V 

YRk = ^ Qn, 

1=1 


(k = N B \,N B i) 


The model for the auxiliary bearings is shown in Figure 2. The governing equations of motion are: 


A/efcA" Bk + C'skY Bl; + A Bk-Y Bk 

= F nk cosa k + \l Bk g, ( 2 - a ) 

M Bk^Bk + C Bk Yg k + K BkY Bk 
= F nk 6in a k , 


( 2 . 6 ) 
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where 


ak — tan 


-1 Ynfe - Y) 


Bk 


A Rk ~ A ' Bk 
(k = A'bi, iV/30) 


Ay A: = - fnit COS C** + A/fc g + A/*. e fi 2 COS (fif), 

3 ~ Fnk sin ock + Mk e fi 2 sin (fit). 

The rotor/bearing interaction is represented with the normal force F n k 

r _ / Kc^k* A* > 0, 

n *~l0, At < 0, 

where 

A* = (X/?* - -Y Dk)cosa k -f (Yji k - } r Bk )sina k - 6 k . 


( 3 ) 


DISCUSSION OF RESULTS 

The model described in the previous section was used to study the steady-state dynamic behavior 
of such a system using direct numerical integration of the equations and a harmonic balance code. 
Studies were performed for varied values of imbalance, support stiffness, and support damping, 
respectively. For the purposes of the current work, the two auxiliary bearings are located at nodes 
3 and 33, respectively. It is assumed that they are identical in terms of stiffness, damping and 
friction characteristics. The nominal system parameters used for the simulation study are M bi =0.0023, 
M b 2=0.0024, K C = 2.855e+6, and £=0.03. All the results that are presented correspond to node 3, 
the location of bearing 1. 

The dynamical behavior of a rotor supported by bearings with clearance coupled with the nonsymmetry 
resulting from gravitational effects can be quite complex. The harmonic balance method is first used 
to investigate the synchronous behavior for small imbalance values. (Please note that the complex 
frequency contents associated with medium and large imbalance values makes it a formidable task 
to apply the harmonic balance method for other cases.) Figure 3 is a typical plot of the steady 
state response amplitudes as functions of the rotor speed with and without bearing clearance, flic 
nonsymmetric effects resulting from gravity loading is similar to what occurs for a rotor supported by 
nonsymmetric bearing stiffnesses with regard to influence on critical speeds. The first critical speed 
splits into two distinct values. For the A direction, the gravity force tends to keep the rotor in contact 
with the bearing at low operating speed. The apparent stiffness in this direction is approximately the 
same as K q and the resulting critical speed about the same as the critical speed for the linear case 
((5=0). For the Y direction, the clearance results in a lower apparent stiffness and a lower critical speed. 
Several higher order pseudo-critical speeds are also created in the operating speed range (about 1500 
rad. /sec). It should be noted that the response in the A" direction also departs from the linear case 
at high operating speed. This is because the imbalance force becomes dominant at high rotor speed 
which in turn makes the gravity force less significant and the clearance effect more important. 

Figure 4(a) shows some typical results for varying imbalance. Some of the more dynamically interesting 
results occur for cases of relatively large imbalance. It is observed that imbalance may influence the 
frequency content of the rotor responses quite dramatically at certain operating speeds. There exists as 
many as eight ranges of imbalance values that result in eight different types of rotor responses. In fact, 
subharmonic responses from fi/2 through fi/10 are observed. Those subharmonics are not directly 
related to the s\siom s natural frequencies as were the cases with other researchers' findings (such as 
Ehrich, 1988). Moreover, several types of subharmonic responses may occur for identical parametric 
configurations, but different imbalance values. 

Figures 4(b) and 1(c) show typical results using A q or C u as the variable parameter. Clear routes 
to chaos are observed. A- !\ b increases beyond certain value, a period-doubling bifurcation always 
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FIGURE 3 INFLUENCE OF CLEARANCE 
ON CRITICAL SPEEDS 
(e=0.0001, C b = 150.0. 
K b =0.313x10 6 ) 


FIGURE 4(a) BIFURCATION DIAGRAM 
WITH IMBALANCE AS THE 
VARIABLE PARAMETER 
(12=1000, <5=0.002, C b = 150. 
K b =0.313x10 6 ) 


takes place. As C g decreases below certain value, the responses always become chaotic. However, the 
bifurcation type is not well defined as for varied Kg. It should be pointed out that even though a lower 
Kg may leads to a better system response, it may also fail to protect the magnetic bearings due to the 
fact that it could result in a larger rotor orbit-center offset. In addition, higher bearing damping does 
not necessarily result in synchronous responses. It is observed that, with small imbalance (e < 0.0005), 
the responses are always synchronous for all speeds in the speed range considered (< 1800) if Cg is 
large enough (C B =250.0) and the clearance is small enough (6 < 0.001). As with the other cases, 
higher bearing damping or lower bearing stiffness or lower imbalance tends to increase the probability 
of synchronous responses. 
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FIGURE 4(b) BIFURCATION DIAGRAM 

WITH BEARING STIFFNESS 

AS THE VARIABLE 

PARAMETER 

(11=1500, <5=0.002, C b = 150, 

e=0.0008) 


FIGURE 4(c) BIFURCATION DIAGRAM 
WITH BEARING DAMPING 
AS THE VARIABLE 
PARAMETER 
(0=1500, 6=0.002, e=0.0015, 
K b =0.313x10 6 ) 
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CONCLUSIONS 

As a summary of the results discussed above, the following conclusions can be drawn: 

1. Imbalance may serve to dramatically alter the frequency contents of the rotor responses at certain 
operating speeds. This is particularly evident for cases of large clearance, high bearing stiffness and 
low bearing damping. 

2. For sufficiently high imbalance: 

i. There is a threshold level of damping below which complex dynamical behavior can be expected. 

ii. There is a threshold level of stiffness above which complex dynamical behavior can be expected. 

3. Clear routes to chaos are observed. As bearing stiffness increases beyond a certain value, a period- 
doubling bifurcation takes place which leads to chaos. As bearing damping decreases below a certain 
value, chaos also tends to occur but the bifurcation type is not so clearly defined. 
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